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ABSTRACT 

We present a multi-wavelength photometric study of ~ 15,000 resolved stars in the nearby spiral galaxy M83 
(NGC 5236, D = 4.61 Mpc) based on Hubble Space Telescope Wide Field Camera 3 observations using four 
filters: F336W, F438W, F555W, and F814W. We select 50 regions (an average size of 260 pc by 280 pc) in the 
spiral arm and inter-arm areas of M83 and determine the age distribution of the luminous stellar populations in 
each region. This is accomplished by correcting for extinction toward each individual star by comparing its colors 
with predictions from stellar isochrones. We compare the resulting luminosity- weighted mean ages of the luminous 
stars in the 50 regions with those determined from several independent methods, including the number ratio of 
red-to-blue supergiants, morphological appearance of the regions, surface brightness fluctuations, and the ages 
of clusters in the regions. We find reasonably good agreement between these methods. We also find that young 
stars are much more likely to be found in concentrated aggregates along spiral arms, while older stars are more 
dispersed. These results are consistent with the scenario that star formation is associated with the spiral arms, and 
stars form primarily in star clusters and then disperse on short timescales to form the field population. The locations 
of Wolf-Rayet stars are found to correlate with the positions of many of the youngest regions, providing additional 
support for our ability to accurately estimate ages. We address the effects of spatial resolution on the measured 
colors, magnitudes, and age estimates. While individual stars can occasionally show measurable differences in the 
colors and magnitudes, the age estimates for entire regions are only slightly affected. 

Key words: galaxies: individual (M83, NGC 5236) - galaxies: stellar content 


1. INTRODUCTION 

Understanding the properties of stars and the history of star 

formation in galaxies remains one of the most fundamental 

subjects in astrophysics. The Hubble Space Telescope (HST) 

provides an important tool for this endeavor, since it enables the 

detailed study of stars and star clusters, not only in the Milky 

Way and its nearest neighbors, but in galaxies well beyond 

the Local Group. The Wide Field Camera 3 (WFC3) provides 

a particularly valuable new panchromatic imaging capability, 

with spectral coverage from the near-UV to the near-IR. This 

capability is especially useful for studying the stellar populations 

in nearby galaxies where individual stars are resolved. 


A good example of a multi-wavelength survey of nearby 
galaxies that employs resolved stars and star clusters is the 
ACS Nearby Galaxy Treasury (ANGST) program (Dalcanton 
et al. 2009). This study includes ~65 galaxies out to ~3.5 Mpc 
and provides uniform multi-color (BVI) catalogs of tens of 
thousands of individual stars in each galaxy. Although the 
ANGST program provides an excellent survey for a wide range 
of studies, it does not provide observations in the U band since 
the observations were obtained before WFC3 was installed 
on HST. The U band is particularly useful for age-dating 
populations of young stars, which is the focus of the current 
paper. PHAT (Panchromatic Hubble Andromeda Treasury) is 
a related study that does take advantage of the new (/-band 
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capability on WFC3. While nowhere as extensive as the Multiple 
Cycle Treasury Program PHAT survey, the current study of 
M83 is complementary in the sense that it provides similar 
observations of the resolved stellar component of a nearby 
spiral galaxy and uses quite different analysis techniques, as 
will be discussed in Sections 2 and 3. Future observations of 
M83 (ID: 125 13, PI: William Blair) will expand the WFC3 data 
set available for M83 from 2 to 7 fields. 

M83 (NGC 5236), also known as the “Southern Pinwheel” 
galaxy, is a slightly barred spiral galaxy with a starbursting 
nucleus located at a distance of 4.61 Mpc, i.e., (m — M ) o = 
28.32 ± 0.13 (Saha et al. 2006). The Ha emission can be used 
to pinpoint regions of recently formed massive stars along the 
spiral arms, while red supergiants can be found throughout the 
galaxy. The color-magnitude diagram (CMD) of resolved stars 
is a powerful diagnostic tool for understanding the stellar evo- 
lution and history of star formation of galaxies in detail. By 
comparing stellar evolution models to observed CMDs, we are 
able to determine the ages of the stellar populations in galaxies. 
However, in active star-forming regions, the stars are often par- 
tially obscured by dust. Applying a single extinction correction 
value for the whole galaxy often results in overestimates or un- 
derestimates of the ages of individual stars. If the CMD is the 
only tool used to determine the ages, the spatial variations of 
dust extinction in galaxies and their effects on the determined 
ages are not readily apparent. The additional information avail- 
able from the color-color diagram can remedy this problem. By 
using techniques developed in this paper, we can constrain the 
variations of dust extinction across M83 and make corrections 
for individual stars. We also focus on spatial variations in the 
stellar properties throughout the galaxy. 

The HST WFC3 observations of M83 were performed in 
2009 August as part of the WFC3 Science Oversight Committee 
Early Release Science (ERS) program (ID: 11360, PI: Robert 
O’Connell). The central region (3.2 kpc x 3.2 kpc) of M83 was 
observed in 2009 August. A second adjacent field to the NNW 
was observed in 2010 March and will be included in future 
publications. Details of the WFC3 ERS data calibration and 
processing are given by Chandar et al. (2010) and Windhorst 
et al. (2011). 

In this paper, we present UBVI photometry of resolved stars 
in M83 and the resulting color-color and color-magnitude 
diagrams. We use (F336W — F438W) versus (F555W — 
FSIAW) color-color diagrams to constrain the variation of dust 
extinction along the line of sight of each individual star, and 
we correct F814W versus (F555Vk — F814W) CMDs for the 
extinction of individual stars to determine their ages. The high 
sensitivity and the superb resolving power in the WFC3 F336W 
band play a key role in allowing us to develop our extinction 
correction techniques and demonstrate the performance of the 
newly installed WFC3. This results in improved stellar age 
estimates and a better understanding of the recent star formation 
history of M83. 

Our investigations are focused on the following: 

1 . Do we see spatial variations of stellar ages in M83 ? If so, can 
we use these variations to learn more about the evolution 
of the galaxy and what triggers star formation? 

2. How well do various age estimates correlate (i.e., resolved 
stars, integrated light, clusters, number ratio of red to blue 
stars, Ha morphology, stellar surface brightness fluctua- 
tions, presence of Wolf-Rayet stars)? 

This paper is organized as follows. In Section 2, we describe 
the observations, photometric analysis, and the extinction cor- 


rection method. CMDs and color-color diagrams of the 50 se- 
lected regions are then used to provide stellar age estimates, 
as described in Section 3. Comparisons with other age esti- 
mates (i.e., integrated light and star clusters) are also made. In 
Section 4, we compare the stellar ages to a variety of other pa- 
rameters that correlate with age, including red-to-blue star ratios, 
Ha morphology, and surface brightness fluctuations. Section 5 
includes a discussion of how these comparisons can be used as 
diagnostics, with special attention paid to the question of what 
spatial variations can tell us about the star formation history in 
M83, and on the question of how star clusters might dissolve 
and populate the held. A summary of our primary results is pro- 
vided in Section 6. In Appendix A, we investigate the locations 
of sources of Wolf-Rayet stars in our M83 held and discuss the 
correlation with the positions of young regions. In Appendix B, 
we describe how the spatial resolution affects our measured col- 
ors and magnitudes and our age estimates of stars in our M83 
held. 

2. HST WFC3 OBSERVATIONS 
2.1. The WFC3/UVIS Data 

The data were obtained in seven broadband and seven 
narrowband hlters in the WFC3 UVIS and IR channels. Details 
are described in Dopita et al. (2010) and Chandar et al. 
(2010). In the current paper, we use four broadband images. 
The hlters and exposure times are F336W (1890 s), F438W 
(1880 s), F555W (1203 s), and F814W (1203 s). Although 
we do not make transformations to the Johnson-Kron-Cousins 
UBVI system, we will refer to these hlters as U. B, V, and 
I for convenience. Three exposures were taken in each hlter 
at different dithered positions to remove cosmic rays, the 
intrachip gap, and to partially compensate for the undersampled 
point-spread function (PSF) of the WFC3 /UVIS channel. The 
raw data were processed using the MULTIDRIZZLE software 
(Koekemoer et al. 2002) with an effective pixel size of 0'.'0396, 
which corresponds to 0.885 pc pixel 1 at the distance of 
4.61 Mpc (Saha et al. 2006). The resulting multidrizzled images 
are combined, cosmic-ray removed, and distortion-corrected. A 
color composite of the F336W, F555W, F814W (broad bands) 
and F502N, F657N (narrow bands) images is shown in Figure 1 . 
It covers the nuclear region of M83, part of its eastern spiral 
arm, and the inter- arm region. Stars in the cores of compact 
star clusters are not resolved at this resolution, but the brighter 
stars in the outskirts of the star clusters and in the field are 
generally resolved into individual stars. Details of the effects 
of spatial resolution are discussed in Appendix B. Since we 
are interested in spatial variations of stellar ages in M83, we 
selected 50 regions in the spiral arm and inter-arm areas. Boxes 
outlined in blue (very young — age< 10 Myr, see Section 3), 
yellow (young — 10 < age <20 Myr), and red (intermediate- 
aged — age > 20 Myr) in Figure 1 show the 50 regions selected 
for detailed study in this paper. 

2.2. Photometry and Artificial Star Tests 

Photometric analysis of the WFC3 M83 data was performed 
on the U, B, V , and I images using the DoPHOT package 
(Schechter et al. 1993) with modifications made by A. Saha. 
Additional routines to derive parameters for the kurtosis of the 
analytic PSFs used by DoPHOT, and for post-processing DoPHOT 
output to obtain calibrated aperture-corrected photometry, were 
performed using customized IDL code written by Saha. A more 
detailed description of these procedures can be found in Saha 
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Figure 1. Color composite of the M83 WFC3 images [Image Credit: Zolt Levay (STScI)]. The F336W image is shown in blue, the F502N ([O hi]) image in cyan, the 
F555W image in green, and the combined F814W and F657N (Ha) image in red. The 50 selected regions are outlined in boxes of three different colors based on the 
values of Age CMDa uto determined in this paper (see Section 3.1 for details). The regions with ages in the range 1-10 Myr are outlined in blue, ages of 10-20 Myr in 
yellow, and ages greater than 20 Myr in red. 


et al. (2010). The kurtosis terms were derived by fitting the 
functional form of the DoPHOT analytic PSF to 20 relatively 
bright and isolated stars. The remaining shape parameters (cr x , 
<jy, and o X y) are dynamically optimized within DoPHOT. 

The difference between the shape parameters for an individual 
object and those for a typical star was calculated and used to 
classify the object as a star, a galaxy, or a double star. The 
shape parameters were determined separately for each image. 
The numbers of objects classified as stars are approximately 
20,000 in U- and B-, 17,000 in V-, and 15,000 in /-band 
images, respectively. The calculated FWHMs are U — O'.' 096 
(2.42 pixel), B = 0'.'093 (2.34 pixel), F= 0T091 (2.29 pixel), and 
/ = 0? 1 1 3 (2.85 pixel). All magnitudes are in the WFC3-UV1S 
VEGAMAG magnitude system, calculated using Equation (4) 


in Sirianni et al. (2005) and the latest zero-point magnitudes: 
U — 23.46, B — 24.98, V = 25.8 1 , and I = 24.67 mag, provided 
by STScI at the HST/ WFC3 Web site. 22 

A multi- wavelength UBVI catalog was then constructed by 
matching the individual objects in each catalog between these 
filters. 

Interestingly, we found that the extreme image crowding 
in these images, combined with the excellent panchromatic 
sensitivity and unprecedented spatial resolution provided by 
the newly installed WFC3, introduced a new challenge during 
this matching step, with large numbers of “pseudo” matches 


22 http://www.stsci.edu/hst/wfc3/ 
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Figure 2. Image cutouts (1572 x 1379, 340 pc x 310 pc) of Regions 26 and 27 in the F555W (top) and the color-composite images (bottom). Stars circled in red 
have (V — 7)>0.8 mag. The “postage stamp” images on the left are “pseudo” matched candidates, while the images on the right appear to be heavily reddened (see 
Section 2.2 for details). 


occurring if too large a match radius was used. We discuss the 
effects of spatial resolution in Appendix B. 

Figure 2 demonstrates this problem by showing the F555W 
and color composite image cutouts near Regions 26 and 27. 
In the middle panel of Figure 2, stars circled in red are objects 
with (V — /) colors that are very red ((V — /) > 1.2 mag) but 
have ( U — B ) values that are fairly blue (—0.4 < (U — B) < 
— 1.7 mag). These can be seen as the spray of points in the upper 
right parts of the top right panels of Figures 3 and 4, which will 
be discussed below. As shown by the postage stamp images in 
Figure 2, about half of these are likely to be blue stars with very 
high values of reddening (i.e., the bottom panel and the small 
postage stamps on the right show that they are in dusty regions), 
while the other half are very close superpositions of at least one 
red and one blue star, which we will call “pseudo” matches (i.e., 
the panels on the left show a strong blue and red gradient across 
many of the objects). 

We found that using a very small matching radius of 0.5 pixels 
was needed to minimize “pseudo” matches in our catalog. 


Earlier attempts using a matching radius of 3 pixels resulted 
in ~10% pseudo matches. Accurate matching also requires 
precise geometric distortion corrections in all filters. However, 
even when these conditions are met, there is a finite chance that 
two different stars will fall within the same aperture. This is 
most easily seen by blinking the U- and /-band images. While 
the single stars stay in the same position, the pseudo match stars 
move slightly, showing they are two different stars with different 
colors. 

The final stellar catalog contains ~ 15,000 objects. We find 
that our procedure results in only a few percent of pseudo 
matches, as will be discussed in Section 2.4 (i.e., the upper 
right panel of Figure 4). 

To measure the photometric completeness, we performed 
artificial star tests with the same detection and photometry 
procedure as applied to the actual stars. We inserted 1000 
artificial stars with Gaussian PSFs and FWHMs appropriate 
for these stars at random positions into all four images. The 
magnitude of the inserted stars varied from 22 to 28 mag 
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Figure 3. Top: the color-magnitude diagram (CMD) and color-color diagram of all stars in our M83 image corrected for Galactic foreground extinction but not for 
internal extinction. Padova isochrones of logr (age/yr) = 6.55, 6.85, 7.0, 7.2, 7.5, 7.8, and 8.0 for a metallicity of Z = 0.03 (1.5 Zq) are overlaid in both panels. 
The dashed line in cyan represents the 50% photometric completeness level. The arrow in each panel indicates the direction of the reddening vector with Ry = 3.1. 
Magnitudes are on the Vega scale. The distance modulus (m — M ) o = 28.32 mag is used to calculate the absolute magnitude Mpsuw- Bottom: histograms of the 
distribution of ages and masses, as determined from the CMD (see Section 3.1 for details). The age quoted in the bottom left panel is the luminosity-weighted mean 
age. 


in steps of 0.25 mag. The 50% completeness levels at 5 a 
detection thresholds for a typical region are reached at recovery 
magnitudes of approximately 24.3, 25.1, 25.3, and 24.8 mag 
for the U, B, V, and I images, respectively. These limits can 
vary by about a magnitude, depending on the brightness of the 
background in a given region. The exception is Region 48, which 
includes the nucleus. The completeness thresholds are two or 
more magnitudes brighter in this region due to the very high 
background. This region has therefore been excluded from the 
rest of the discussion, but it is included in Table 1 for reference. 

2.3. CMDs and Color-Color Diagrams in M83 

The top panels of Figure 3 show the I versus (V—/) CMD 
(left) and the (U — B) versus (V — I) color-color diagram (right) 


of resolved stars extracted from the cross-matched catalog. 
The CMD shows the presence of young main-sequence (MS) 
stars, transition or He-burning “blue loop” stars, red giant stars 
(hydrogen shell-burning and hydrogen + helium shell-burning 
stars), and a relatively small number of low-mass red giant 
branch (RGB) stars. This is because the tip of the low-mass 
RGB stars is at I = 24.7 mag (Karachentsev et al. 2007), which 
is roughly the same as our completeness threshold. The ages of 
these stars range from s»l to ~ 1 00 Myr. The Padova isochrones 
(Marigo et al. 2008) ranging in age from log r (age/yr) = 6.55 
(3.5 Myr) to log r (age/yr) = 8.0 (100 Myr) are included in 
Figure 3 for reference. We note that the younger isochrones 
(e.g., 1 Myr) fall nearly on top of the 3.5 Myr isochrone and 
hence are not included in the diagram. 
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Figure 4. Top: the CMD and color-color diagram of all stars corrected for the internal extinction determined for each individual star. The blue and two red dotted 
lines indicate boundaries for the locations of stars in the color-color diagram uncorrected for internal extinction (see Section 2.4 for discussion). Two arrows in the 
color-color diagram show reddening vectors: one with the standard reddening vector {Ry = 3.1) and the other with the flatter reddening vector (Ry = 5.7), which 
appears to be more appropriate for some of the data in the green triangle (see Section 2.4 for discussion). The dashed line in magenta shows the Mj < —5.5 mag 
cutoff used for the Age CMDauto estimate. Bottom: same as described in Figure 3. 


The arrows in both panels show Ay = 1 reddening vectors 
for M83 using the Ry = 3.1 extinction curve of Cardelli et al. 
(1989). Corrections have been made for foreground Galactic 
extinction (Schlegel et al. 1998) using A F i 36 w = 0.361, 
Af438iv = 0.290, Af555w = 0.229, and Af 8 i 4 iy = 0.133 mag, 
respectively. 2 ’ In Figure 3, we only plot stars with photometric 
errors less than 0.25 mag in the filters used for each diagram. 
The numbers of objects plotted in the CMD and the color-color 
diagram are about 12,000 and 8500 stars, respectively. Since 
approximately 30% of the stars are not detected in one or more 
of the U, B, or I filters, due to either reddening or intrinsically 
red or blue colors, fewer stars are plotted in the color-color 
diagram. 

We adopt a value of 1.5 times solar metallicity (Z = 0.03), 
the highest metallicity isochrones available from the Padova 

23 NASA/IPAC Extragalactic Database (NED): http://ned.ipac.caltech.edu/ 


database, for M83 based on Bresolin & Kennicutt (2002). A 
test using solar metallicity (Z = 0.019) isochrones showed 
that the ages for a typical region would be «s2 Myr older if 
the lower metallicity is adopted. This agrees with the study 
by Larsen et al. (2011), who found that the simulated CMDs 
of young massive clusters in M83 with solar and super- 
solar metallicity isochrones would not look much different. 
Theoretical isochrones calculated for the WFC3 filters and Z = 
0.03 (1.5 Z Q ) from the Padova database 24 (Marigo et al. 2008) 
are overlaid on both panels in Figure 3. The dashed line in cyan 
shows the 50% completeness level in I and ( V — I), using the 
completeness threshold numbers from the artificial star test in 
Section 2.2. 

The bottom panels of Figures 3 and 4 show the histograms 
of ages and masses for the same stars plotted in the CMD 

24 CMD version 2.3; http://stev.oapd.inaf.it/cgi-bin/cmd 
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Region 48 is the nucleus of M83. It is not included in the analysis since the bright background results in a completeness limit that is more than 2 mag brighter than other regions. 
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and color-color diagram. The ages and masses of stars were 
estimated by finding the closest match between the / and (V — I) 
values for each star, using a fine mesh of the stellar isochrones 
plotted in the CMD and color-color diagram in Figures 3 and 4. 
More details about the stellar age-dating from the CMDs and 
the stellar isochrones are given in Section 3.1. 

2.4. Extinction Corrections 

Since M83 has an intricate structure of dust lanes associated 
with active star-forming regions in the spiral arms and thin 
layers of dust in the inter-arm area in Figure 1, we cannot use 
a single value of internal extinction and apply it to correct for 
the extinction of all stars in a given region. In the color-color 
diagram in Figure 3, we notice that while the bluest stars match 
the model isochrones quite well after we correct for Galactic 
foreground extinction, the majority of stars are found redward 
of the models, implying a typical extinction of Ay ~ 0.5 mag 
(for the densest part of the data “swarm”). As described below, 
we can determine the reddening of these stars if we assume 
that they are intrinsically blue stars that belong on the Padova 
models (i.e., with (V — I) — —0.3 mag). A visual inspection 
supports this interpretation. Stars with observed values of 
(V — /) —0.3 mag are found in areas with no obvious dust 

features surrounding them, while stars with (V — I) > 0.5 mag 
are near dust filaments. 

Our basic approach to correct the colors (and luminosities) 
of individual stars for the effects of extinction is to backtrack 
the position along the reddening vector for each data point, until 
we hit the stellar isochrones in the color-color diagram. What 
makes this method work is the fact that all the isochrones are 
nearly on top of each other in the range (U — B) < 0.0 mag. 
The way we apply this in practice is to match the observed and 
predicted values of the reddening-free Q parameter defined by 

E(U - B) 

2ubvi = (U — B) — -±- -J x (V - /) (1) 

E(V - I) 

(Mihalas & Binney 1981). Using observed values of UBVI, we 
compute the observed Qubvi values by using the standard slope 
of E(U - B)/E(V - I) = 0.58 (Mihalas & Binney 1981; 
Whitmore et al. 1999), which corresponds to Ry = 3.1. The 
predicted Q values are calculated using the 1.5 Z Q Padova 
isochrones. A complication is that we cannot simply use the 
matched Q values to determine corrected values of ( U — B ) 
and ( V — /), since they would, by definition, fall precisely 
on the isochrones (i.e., we would be using circular reasoning). 
Instead, we solve for the extinction in ill — /), which is partly 
independent and uses the longest wavelength baseline. We then 
calculate the extinction values in UBVI and the color excess 
values in (U — B) and (V — I). The averages of the corrected 
internal extinction and color excess values by our star-by-star 
correction method are Ap^w = 0.696, A F aizw = 0.559, 
^F 555 iy = 0.441, A/r 8 i 4 w = 0.256 mag, and E(U — B ) = 0.137 
and E(V — /) = 0.185 mag. This is similar to the average 
internal extinction (Ap^w = 0.425 mag) estimated from the 
cluster spectral energy distribution (SED) fitting by Bastian 
et al. (2011). However, our results differ from the average 
extinction {Ay = 0.671 mag) of 45 clusters in the nucleus 
(~20" in diameter) of M83, determined from Ho'/H/J ratios 
by Harris et al. (2001). This is consistent with the finding of 
larger extinction in the nuclear region by Chandar et al. (2010). 

In the color-color diagram (top right panel) of Figure 3, we 
note that while most of the data points are consistent with 
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the standard reddening vector, the stars located in the green 
triangle appear to follow a flatter reddening vector. Whether 
this is actually due to a different reddening law in M83 (e.g., 
as has been suggested for heavily extincted regions such as 4, 
12, and 20) or is some sort of artifact (e.g., due to photometric 
uncertainties, unresolved star clusters, or the “pseudo” matching 
problem discussed above) is difficult to determine from our 
present data set. However, for our specific needs the precise 
answer to this question is not critical. Our approach will be to 
correct the data points in the triangle back to a position near the 
top of the isochrones, as shown by the intersection of the red 
dotted line and the isochrones in Figure 4. This is equivalent to 
using a range in Ry between 3.1 and 5.7 (the value represented 
by the red dotted line). 

The upper panels of Figure 4 show the CMD and color-color 
diagram corrected for the internal extinction in M83. The arrows 
in both panels are the reddening vectors. We do not apply an 
extinction correction for stars with (V — I) ^ 1.2 mag or above 
an extrapolation of the flatter extinction vector from the bluest 
possible isochrones (see the red dotted lines in the upper right 
panel of Figure 4), since we believe that the colors of many 
of these sources are incorrect (i.e., roughly half are likely to 
be pseudo matches, as discussed in Section 2.2 and shown in 
Figure 2). We note that only about 4% of the stars fall in this part 
of the diagram. These stars are removed from the subsequent 
analysis. If we include these stars, our age estimates typically 
increase by about 1 Myr. 

Also, stars below the blue dotted line in Figure 4 are not 
corrected for extinction, since these are MS turnoff stars (i.e., 
blue loop or transition stars in general), and hence do not have 
unique UBVI colors for the stars along a given reddening line, 
unlike the MS stars with (U — B) < 0.0 mag. 

Based on the extinction-corrected CMD in Figure 4, we 
note that (1) the corrected data now show a relatively narrow 
distribution of points along the left side of the CMD in agreement 
with ages of ~3 Myr and (2) while the data now align with the 
model isochrones much better, there is still a fair amount of 
scatter (especially for I < 24), primarily due to observational 
uncertainties from the fainter stars. 

The bottom panels of Figure 4 show the distribution of 
ages and masses of stars after we apply our extinction cor- 
rection technique. Based on the automatic method described in 
Section 3. 1, we note that the luminosity-weighted mean ages de- 
termined from the corrected colors and luminosities are signifi- 
cantly younger (~14 Myr; Figure 4) than the mean ages deter- 
mined from the uncorrected colors and luminosities (~27 Myr; 
Figure 3), as we expected. 

3. AGE-DATING POPULATIONS IN M83 

In this section, we discuss several different methods of age- 
dating stellar populations in M83, using both individual stars 
and star clusters, and then intercompare the results. We do not 
necessarily expect field star and cluster ages to agree, since the 
dissolution of clusters may bias the observed cluster population 
toward younger ages than the surrounding field stars. 

3.1. Age-dating the Resolved Stellar Population Using 
Color-Magnitude Diagrams 

As shown in Figure 1, we selected 50 regions that cover 
spiral arm and inter-arm areas in order to study the recent star 
formation history of M83. Figures 5 and 6 show cutouts of 
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Figure 5. CMDs, color-color diagrams, and the image cutouts of Regions 12 and 7. Data points before and after the individual extinction correction are plotted as 
black plus signs (upper: uncorrected) and blue crosses (lower: corrected). The Padova isochrones in these diagrams are the same as the ones plotted in Figure 3. See 
Section 2.4 for details. 
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F555W - F814W 




Figure 6. Same as Figure 5 for Regions 29 and 2. 
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log (Age HK [yr]) 



log (Age c 


„[yr]) 


Figure 7. Correlation plots between the ages of the stellar populations in the 50 regions from CMD isochrone fitting. Left: correlation between the manual ages 
determined by HK and BCW. Right: correlation between the average of manual ages and automatic ages (see Section 3.1 for details). The best linear fits are shown in 
black with the significance of the correlation in units of o (i.e., slope/uncertainty) in the top left of each panel. The dashed line is the unity line. 


Regions 12, 7, 29, and 2. 15 These figures include color-color 
diagrams of the stars in these regions, as well as CMDs that are 
uncorrected (upper) and corrected (lower) for extinction. These 
four regions highlight the range in age of the dominant stellar 
population, from very young to intermediate ages, i.e.. Region 
12 with very strong Ha emission superposed on the cluster stars; 
Region 7 with a large bubble of Ha emission surrounding the 
stars and clusters; Region 29 with no Ha emission, but large 
numbers of bright red and blue stars; and Region 2 with no Ha 
emission and fainter stars. 

A cursory glance at the CMDs of these four regions (Figures 5 
and 6) shows clear differences. The primary differences are, 
as predicted by the isochrones, that younger regions contain 
(1) bluer MS stars, (2) larger numbers of upper MS stars, and 
(3) larger ratios of blue to red stars. Our regions do not, however, 
contain only stars of a single age. Even these four regions, which 
were chosen to include a single dominant stellar population, 
contain a mix of young and old stars. In this paper, we are 
primarily interested in the ages of the bright young stars in these 
regions, which dominate CMDs in luminosity-limited samples. 

In practice, this focus on the younger population is carried 
out by imposing a magnitude cutoff of M/ = —5.5 mag (i.e., 
I = 22.82 mag, which corresponds to the age cutoff of ~60 Myr 
in the CMD) for the stars used to estimate ages (i.e., the magenta 
dashed line in Figure 4). This luminosity limit allows us to stay 
above the completeness limit for all but the very reddest stars 
(see Figure 3), while also focusing on evolved stars and those 
on the upper MS, which are most sensitive to the age of a stellar 
population. Our primary goal is to obtain reliable “relative” 
estimates for the 50 regions (rather than absolute ages), and our 
magnitude limit is sufficiently deep to accomplish this goal. 

We compare the extinction-corrected colors and magnitudes 
of stars with the Padova stellar isochrones in two ways to esti- 
mate the age of the dominant stellar population in each region. 
First, two of the authors (H.K. and B.C.W.) independently es- 
timated ages based on a visual comparison of the CMDs and 
isochrones. We focused on features such as the color of the MS 


25 The color-magnitude diagrams and color-color diagrams of all 50 regions 
(similar to Figures 5 and 6) are available on the WFC3 Early Release Science 
Archive Web site (http://archive.stsci.edu/prepds/wfc3ers/). 


stars, the number of stars in the upper MS, and the number ra- 
tio of blue to red stars. The left panel of Figure 7 shows that 
the independent, visually determined age estimates are in good 
agreement with a 7.9a (slope/uncertainty of the best linear fit) 
correlation and a slope within la of the unity value. The average 
ages of these two independent (manual) estimates are listed in 
column Age CMDman in Table 1 . 

Next, we estimated the ages automatically, finding the closest 
match in both age and mass for each star by comparing the 
extinction-corrected I magnitude and (V — I) color with a fine 
mesh of stellar isochrones generated from the Padova models 
(logr (age/yr) ranging from 6.05 to 8.35 in steps of 0.05). We 
then calculate a luminosity-weighted age from the distribution of 
individual stellar ages in each region. The result of the automatic 
age estimate is listed in column Age CMDaut0 in Table 1 . 

To assess the uncertainties in our age and mass determinations 
for each star, we ran a test for four different cases by adding 
and subtracting the photometric errors in the (V — I) color and 
/-band magnitude. This moves each object on the CMD to the 
right (+( V - /) error X left (— (V - /) er ror), down (+m /cn J, and up 
(— m / error ). We then estimated the ages for these four cases in the 
same way as we did to get Age CMDaut0 . The results of this test 
are summarized in column AAge CMDauto in Table 1. As shown 
in the AAge CMDauto column in Table 1 , the variation in the age 
estimates is typically only a few Myr (i.e., <10%). 

Figure 8 shows the age and mass histograms of individual 
stars in the four regions shown in Figures 5 and 6. Luminosity- 
weighted mean values of the stellar ages are shown in the upper 
right of the left panels. The age sequence of these regions is 
as expected from the youngest in Region 12 to the oldest in 
Region 2. 

One unexpected result was that while Region 2 was selected 
to be representative of an older field population, it also contains 
a fair number of bright young stars. This suggests that young 
individual stars may form in the field in low numbers. Another 
possibility is that many of these are “runaway” stars, i.e., 
young stars that have been dynamically ejected from their birth 
sites within young star clusters due to interactions with other 
stars. Whitmore et al. (2011a) also suggested this as a possible 
explanation for isolated H ii regions ionized by a single, massive 
star, i.e., “Single-Star” Hu (SSHII) regions (see Figure 3 from 
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Figure 8. Histograms of the distribution of ages and masses of the individual stars in Regions 12, 7, 29, and 2, as determined by the automated isochrone fitting 
technique (see Section 3.1). The age indicated at the top right in each age histogram is the luminosity-weighted mean age of all stars with Mj < —5.5 in each region. 
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Figure 9. Image cutout (18" 3 x 16"9, 410 pc x 380 pc) of the area near Region 
9 as an example of a “field” region with old red giant stars and a small population 
of isolated young blue stars. See Section 3.1 for details. 


Whitmore et al. 201 lb). Either explanation would be interesting 
and should be investigated in a future study. Figure 9 shows 
a similar “field” Region 9 with primarily old red giant stars 
(barely discernible), but which also includes a small population 
of isolated young blue stars and four compact clusters ranging 
in age from intermediate to old. 

The right panel of Figure 7 compares the Age CMDtnan and 
Age CM Dauto estimates for the different regions. A unity line is 
included for comparison and shows that there is good agreement 
(7.6 a correlation) between the two methods that varies from 
the unity line by only la. We use the objective luminosity- 
weighted age estimates (Age CMDauto ) as our preferred method 
for estimating ages of the resolved stellar populations in the 
remainder of this paper. 

3.2. Comparison with Previous Age Estimates 

Larsen et al. (2011) recently estimated the ages of two 
star clusters in M83 by comparing photometry of individual 
stars in the outskirts of the clusters with isochrones. These 
clusters are within Regions 3 and 11 here. Larsen et al. (2011) 
used a more sophisticated method, comparing observed and 
predicted magnitude distributions for red and blue stars using 
Hess diagrams. They did not, however, correct the magnitudes 
or colors of their stars for the effects of spatial variation or 
extinction, but since these are intermediate-age clusters with 
minimal internal extinction, this is not a large effect in these 
particular cases. They find log r of 7.44 (27.5 Myr) and 7.45 
(28.2 Myr) for the clusters in Regions 3 and 11. We find 
similar ages of 21.9 Myr (Region 3) and 26.7 Myr (Region 11) 
for the dominant luminous stellar populations, despite the fact 
that our photometry includes a significantly larger area around 
each cluster. Hence, we find good agreement between our age 
estimates and those of Larsen et al. (201 1). 

Each method has its own strengths and weaknesses. Our 
approach has the advantage of correcting the photometry of 
individual stars for the effects of extinction to improve our age 
estimates, but the disadvantage that this requires imaging in 
more than three optical bands with a crucial need for the t/-band 
observation. The Larsen et al. (2011) method does not correct 
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for the effects of extinction but allows for a better determination 
of the star formation history over a wider age range (i.e., out to 
~1 Gyr). 

3.3. Comparison between Age CMDaut0 and 
Compact Cluster Ages 

Much of our recent work has involved age-dating star clusters 
using the integrated light ( UBVI and Ha emission) from the 
cluster and model SED fitting. To compare the age of star 
clusters in each region to the age determined from individual 
stars in each region, we adopt the ages and luminosities of 
compact clusters identified in our previous work (Chandar et al. 
2010). Details of the cluster age-dating method can be found in 
that paper. 

We note that the regions sampled in this study are large (i.e., 
several hundred pc 2 ) and may include stars and star clusters 
spanning a wide range of ages. This requires making the 
comparison between the ages of field stars and star clusters 
with caution. Only the largest star-forming complexes within 
spiral galaxies have such dimensions (e.g., the giant Hn region 
NGC 604 in M33; Freedman et al. 2001). More typical regions 
are much smaller and hence would not dominate the entire 
field. This will tend to weaken our correlations, especially when 
comparing resolved stellar ages with cluster ages. 

Using the ages of star clusters determined from the SED 
fitting, we compute the luminosity-weighted mean age of each 
region, as we did for stars in the same region. The ages of 50 
regions are listed in column Age cl in Table 1 . The upper panel of 
Figure 10 shows a comparison between the luminosity-weighted 
mean age (Age cl ) of the compact clusters in a given region and 
our stellar age (Age CMDauto ) estimates based on the CMDs of all 
the resolved stars in the region (as discussed in Section 3.1). The 
correlation between these two ages is fair (5.4a correlation). We 
note that while the midpoints for the two methods are similar, 
the slope is steeper, as discussed in more detail in the next 
section. 

3.4. Comparison between Age CMDauto and Age Estimates 
from Integrated Light from the Entire Regions 

We obtain an independent estimate of the age of each region 
by measuring the colors of the entire region and perform a simple 
SED fit, in the same way as we did for star clusters (i.e., Chandar 
et al. 2010). The lower panel of Figure 10 shows a comparison 
between our age estimates for the resolved stellar populations 
using CMDs (Age CMDaut0 ) and those from the integrated light 
(Age R ) in each region. 

While there is a fair correlation (4.6a) between the CMD and 
integrated photometric age estimates, there is also a fair amount 
of scatter and an apparent gap in the region ages in the range 
6.8 < log (age/yr) <7.2. This gap is similar to the well-known 
artifact for cluster age estimates, which is due to the looping 
of predicted cluster colors in this age range (see Chandar et al. 
2010). We also note that the slope of the relation is steeper than 
the unity vector, with integrated light ages ranging to much lower 
ages than the CMD ages. This is similar to the comparison with 
Age cl , as discussed above. This is probably caused by a variety 
of effects, including (1) the integrated age estimates take into 
account Ha emission (which is very sensitive to massive young 
stars) while the CMD estimates do not; (2) the 1 Myr isochrones 
are essentially on top of the 3 Myr isochrones, making it rare that 
the minimum age ever gets selected by the software that does 
the matching with the isochrones; and (3) the adoption of an 
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log (Age CMDauto [yr]) 

Figure 10. Top: correlations between the ages of stars (Age CMDa uto) and 
clusters (Age cl ). Bottom: correlation between the ages of stars determined by 
A§ e CMDauto and the integrated light within the region (Age Reg ). See Section 3.4 
for details. The best linear fits are shown in black solid line with 5.4<r (top) and 
4.6<j (bottom) correlations. The dashed line is the unity line. 

My — — 5.5 mag cutoff removes all stars with ages >100 Myr 
from the Age CMDauto estimates, while the light from older stars is 
still included in the integrated light used for the Age Reg estimate 
(hence, the Age CMDaut0 estimate can never be >50 Myr, while 
the Age Reo can be older). 

This difference in slope in Figure 10 highlights the fact that 
each age-dating method has its own idiosyncrasies. By com- 
paring a number of different methods, we begin to understand 
these artifacts better and learn how large the true systematic 
uncertainties can be. 


4. COMPARISON BETWEEN Age CMDaut0 AND OTHER 
PARAMETERS THAT CORRELATE WITH AGE 

We are now ready to compare our age estimates for the 
resolved stellar components in 50 regions to other observables 
that correlate with age. These include (1) number ratio of red to 
blue stars (e.g., Larsen et al. 201 1), (2) morphological categories 
(Whitmore et al. 2011a), and (3) stellar surface brightness 
fluctuations (Whitmore et al. 2011a). 


4.1. Comparison between Age CMDaut0 and Number 
Ratio of Red to Blue Stars 

One property that is expected to correlate with age is the 
number ratio of red to blue stars. At very young age, all stars 
are blue. As the population ages, the number of red giant stars 
(H-shell burning and H+He-shell burning stars) increases. This 
effect can be seen by noting that nearly all of the stars in Region 
12 (Figure 5) are blue, while there are large numbers of red stars 
in Regions 29 and 2 (Figure 6). 

The number ratio of red to blue stars (column red-to-blue 
ratio in Table 1) is calculated by using a criterion that (V — I) 
be redder than 0.8 mag for red stars. The top panel of Figure 1 1 
shows a fairly good correlation (1 .2a) between Age CMDauto and 
the red-to-blue ratio. 

4.2. Comparison between Age CMDauto 
and Morphological Category 

The middle panel of Figure 1 1 shows the correlation be- 
tween our automatic CMD age estimates and the morpholog- 
ical classification for a given region, as defined in Whitmore 
et al. (201 la). Briefly, regions with Ha emission superposed on 
top of the stellar component are type 3 (emerging star clus- 
ters), regions with small Ha bubbles are category 4a (very 
young), regions with large Ha bubbles are category 4b (young), 
and regions with no Ha bubbles are category 5 (intermediate 
age). Based on these criteria, the regions displayed in Figures 5 
and 6 are classified as categories 4a (Region 12), 4b (Region 7), 
5a (Region 29), and 5a/5b (Region 2), respectively. Each region 
was classified independently by two authors (H.K. and B.C.W.). 
The mean of the two determinations is used in what follows and 
is listed in column “Ha Morphology” of Table 1 . 

The middle panel of Figure 1 1 shows that stellar CMD 
age estimates (Age CMDauto ) and morphological categories (Ha 
morphology) show a fair correlation (5.8a), but not as good 
as the strong correlations (9a for log r < 7 and 5 a for log 
r > 7) found with ages for compact clusters in Whitmore et al. 
(201 la). This is probably due to the fact that most of the regions 
are not dominated by a single-age stellar population, making 
the morphological classification for an entire region problematic 
(i.e., there is strong Ha emission in part of the region but none 
in other parts). 

4.3. Comparison between Age CMDaut0 and Stellar 
Surface Brightness Fluctuations 

In Whitmore et al. (2011a), we developed a method to age- 
date star clusters based on an observed relation between pixel- 
to-pixel flux variations (rms) within star clusters and their ages. 
This method relies on the fact that the young clusters with 
bright stars have higher pixel-to-pixel flux variations, while the 
old clusters with smoother appearance have small variations. 
This is because the brightest stars in a 100 Myr population 
have Mi —3 mag, which is below our 50% photometric 
completeness limit (M/ = —3.6 mag) and well below our 
magnitude cutoff of M/ = —5.5 mag used in our stellar age- 
dating method described in Section 3.1. 

The bottom panel of Figure 1 1 shows the correlation between 
the resolved stellar ages (Age CMDauto ) and the surface brightness 
fluctuations (pixel-to-pixel rms) measured in the 50 selected 
regions. We find little or no correlation (1.1a) with a relatively 
large amount of scatter, especially for the younger clusters. Part 
of this scatter may be caused by the two-valued nature of the 
effect as discussed in Whitmore et al. (2011a), with a maximum 
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log (Age CMDaulo [yr]) 

Figure 11. Correlations between the stellar ages (Age CMDauto ) and the number 
ratio of red to blue stars (top), Ha morphology (middle), and pixel-to-pixel 
fluctuations (bottom) for the 50 selected regions. See Section 4 for details. The 
best linear fits are indicated in solid black line with 7 .2a (top), 5.8a (middle), 
and 1.1a (bottom) correlations. 


value of the surface brightness fluctuations around 10 Myr, and 
lower values for both smaller and larger ages. Another reason 
for the scatter is the fact that many of the regions have mixed 
populations, rather than single-age populations. 


5. INSIGHTS INTO THE STAR FORMATION HISTORY IN 
M83 BASED ON SPATIAL VARIATIONS FOR 50 REGIONS 

We first examine how the ages for the 50 regions are 
distributed in Figure 1 . The boxes in this figure are color-coded 
as follows: blue (very young with ages less than about 10 Myr), 
yellow (young with ages between about 10 and 20 Myr), and 
red (intermediate aged with ages greater than about 20 Myr). 
As expected, we find that the very young regions are primarily 
associated with the spiral arms, the young regions tend to be 
“downstream” from the spiral structure (i.e., on the clockwise 
side away from the dust lanes), while the intermediate-aged 
regions are found on both sides of the spiral arms in the inter- 
arm regions. 

A similar, more detailed analysis is being done by 
Chandar and Chien using the cluster ages in M51 and M83 
(R. Chandar & L.-H. Chien 2012, in preparation). They com- 
pare the age gradients they find with models developed by Dobbs 
& Pringle (2010), assuming a variety of different triggering 
mechanisms (i.e., spiral arms, bars, stochasticity, and tidal dis- 
turbances). 

We can perform a similar experiment here, since we have age 
estimates for each individual star. The three panels in Figure 12 
show the distribution of stars with ages in the ranges 1-10, 
15-35, and 40-100 Myr, respectively. The top panel is the 
youngest group of stars, clearly showing that the stars in these 
regions are mostly distributed along the active star-forming 
region (i.e., associated with strong Her emission) in the spiral 
arms. The stars in the middle panel tend to be found slightly 
downstream of the spiral arms, while the older stars are still 
farther out in the inter-arm regions, as expected. 

These diagrams can also be used to examine whether there 
is evidence that most of the young stars form in clusters and 
clustered regions and then dissolve into the field. Pellerin et al. 
(2007) made a similar set of diagrams for the galaxy NGC 1313, 
which supported this interpretation. We find that the young star 
samples show strong clustering, while the older star samples are 
progressively more uniform. Hence, these distribution maps of 
stars in the different age group support the idea that the stars 
form in clusters in spiral arms, and then most of the clusters 
dissolve, populating the field with stars (Lada & Lada 2003). R. 
Chandar & L.-H. Chien (2012, in preparation) will examine this 
subject more quantitatively in the future. 

6. SUMMARY 

Color-magnitude diagrams and color-color diagrams of re- 
solved stars from the multi-band HST /WFC3 ERS observations 
of M83 have been used to measure the ages of stellar popula- 
tions in 50 regions of this well-known face-on spiral galaxy. The 
diagrams show the presence of multiple stellar features, includ- 
ing recently formed MS stars, He-burning blue-loop stars, and 
shell-burning red giant stars. Comparisons between our stellar 
age estimates and a wide variety of other age estimators allow 
us to investigate a number of interesting topics. The primary 
new results from this study are as follows: 

1 . An innovative new technique using a combination of CMD 
and color-color diagrams has been developed in order to 
correct for the extinction toward each individual star and to 
age-date the stellar population. The mean extinction values 
for the 50 regions studied in this galaxy are 0.696 (A f 336 w), 
0.559 (Af438vr)> 0.441 (A^sw), and 0.256 (Afgi 4 W ) mag, 
respectively. 
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Figure 12. Spatial distribution of stars with ages of 1-10, 15-35, and 
40-100 Myr determined in this study. Each panel has the same orientation 
and scale as Figure 1. The top panel is the youngest group of stars, clearly 
showing that the stars in these regions are mostly distributed along the active 
star-forming region (i.e., associated with the strong Ha emission) in the spiral 
arms. The stars in the middle panel tend to be found slightly downstream of the 
spiral arms, while the older stars are still farther out in the inter-arm regions, as 
expected. See Section 5 for details. 
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2. The various age estimators (stellar, integrated light, clusters 
within the region) show fair correlations (i.e., between 4.6a 
and 5.4a). This is expected, since many of the regions have 
a mixture of different-aged populations within them, and 
each technique uses light from a different subset of the 


stars (e.g., the stellar age estimates use only bright stars, 
while integrated light includes all the light). A comparison 
with previous age estimates by Larsen et al. (201 1) shows 
a good agreement (see Section 3.2 for details). 

3. Comparisons between the stellar ages and other parameters 
that are known to correlate with age show a range from very 
good correlations (e.g., with red-to-blue ratios, i.e., 7.2 a) 
and fair correlations (e.g., Ha morphology with 5.8a) to 
little or no correlation (e.g., pixel-pixel rms from Whitmore 
et al. 2011a). This is expected for reasons similar to the ones 
discussed above. 

4. The regions with ages younger than 10 Myr are generally 
located along the active star-forming regions in the spiral 
arm. The intermediate-age stars tend to be found “down- 
stream” (i.e., on the opposite side from the dust lane) of the 
spiral arms, as expected based on density wave models. A 
more detailed comparison with models with various other 
triggering mechanisms (e.g., bars, density waves, tidal dis- 
turbance, and stochasticity; see Dobbs & Pringle 2010) is 
in process by Chandar and Chien for M51 and M83 (R. 
Chandar & L.-H. Chien 2012, in preparation). 

5. As discussed in Appendix A, the locations of Wolf-Rayet 
sources from Hadfield et al. (2005) are in broad agreement 
with the age estimates discussed in the current paper. The 
much better spatial resolution from HST shows that many 
of the Wolf-Rayet “stars” from ground-based observations 
are actually young star clusters. 

6. Effects of spatial resolution on the measured colors, mag- 
nitudes, and the derived ages of stars in our M83 images 
are described in Appendix B. Based on a numerical experi- 
ment using a star cluster NGC 2108 in the Large Magellanic 
Cloud (LMC), we found that while individual stars can oc- 
casionally show measurable differences in the colors and 
magnitudes, the age estimates for entire regions are only 
slightly affected. 

We thank Zolt Levay for making the color images used in 
Figure 1. We also thank the anonymous referee for helpful 
comments. This project is based on Early Release Science 
observations made by the WFC3 Science Oversight Committee. 
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program. Support for program 11360 was provided by NASA 
through a grant from Space Telescope Science Institute, which 
is operated by the Association of Universities for Research in 
Astronomy, Inc., under NASA contract NAS 5-26555. 
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APPENDIX A 

USING WOLF-RAYET STARS TO CHECK FOR 
CONSISTENCY WITH OUR AGE ESTIMATES 

In this paper, we compared the ages resulting from a variety 
of age-dating methods for different star-forming regions in 
M83. Here, we extend our analysis by considering previously 
identified Wolf-Rayet stars in M83, which are ~ 1 -4 Myr old 
(Crowther 2007), as another independent age estimate. 

Hadfield et al. (2005) identified 283 Wolf-Rayet sources in 
M83 from narrowband imaging centered on the Hen 14686 
emission line, plus spectroscopic follow-up. Seventeen of these 
Wolf-Rayet sources fall in our WFC3 field of view. Figure 13 
shows their locations in the HST image, and Table 2 lists the 
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Figure 13. Distribution of Wolf-Rayet sources in our HST WFC3 M83 images. The color image at the center is the UBVIHa composite, and the image cutouts are 
the Ha (left or top) and F814W (right or bottom) images for regions A-E. The top right panel is an original finding chart (Hadfield et al. 2005, A.4684 narrowband 
FORS2 image) for region A. Numbers in blue are the IDs of Wolf-Rayet sources identified by Hadfield et al. (2005). 


Table 2 

Wolf-Rayet Star Candidates in M83 


R.A. 

Decl. 

IDHadfleld 

Reg No. 

A g e Reg( M y r ) 

X 

Y 

Comments 

13:37:00.41 

-29:52:54.1 

72 

- 

- 

3874 

606 

Faint Ha 

13:37:01.18 

-29:52:53.9 

73 

- 

- 

3618 

609 

Faint Ha 

13:37:01.42 

-29:51:25.8 

74 

4 

- 

3512 

2821 

Strong Ha 

13:37:03.06 

-29:50:49.5 

78 

7 

4.2 

2982 

3698 

Faint Ha 

13:37:04.65 

-29:52:47.5 

82 

- 

- 

2710 

772 

Strong Ha 

13:37:04.65 

-29:50:58.4 

86 

12 

2.4 

2477 

3505 

Strong Ha 

13:37:07.35 

-29:49:39.2 

95 

20 

0.8 

1513 

3306 

Strong Ha 

13:37:07.54 

-29:52:53.6 

96 

41 

50.1 

1552 

620 

Strong Ha 

13:37:07.96 

-29:52:07.9 

99 

- 

- 

1418 

1762 

Faint Ha 

13:37:08.25 

-29:51:13.6 

100 

24 

- 

1311 

3116 

Faint Ha 

13:37:08.40 

-29:52:54.9 

102 

- 

- 

1266 

588 

Strong Ha 

13:37:08.53 

-29:52:12.0 

103 

35 

1.1 

1224 

1659 

Strong Ha 

13:37:08.70 

-29:52:28.9 

105 

38 

2.1 

1173 

1239 

Strong Ha 

13:37:08.91 

-29:52:05.7 

106 

31 

4.4 

1101 

1818 

Faint Ha 

13:37:09.02 

-29:52:45.2 

107 

39 

24.0 

1067 

825 

No Ha 

13:37:09.80 

-29:52:36.2 

108 

- 

- 

820 

1049 

Strong Ha 

13:37:10.42 

-29:51:28.0 

109 

26 

4.6 

610 

2756 

Faint Ha 


coordinates of these sources, as well as the region in which they 
are located. 

The lower resolution of the ground-based Hadfield et al. 
(2005) study makes it difficult (in many cases) to uniquely 
identify the Wolf-Rayet sources in the HST image, although we 
can easily determine whether such a source is located within one 


of our 50 regions. Figure 1 3 shows a color image of M83 (center 
panel) with five regions, labeled A through E, that contain all 
17 Wolf-Rayet sources. Here, the Ha and F814W images are 
shown in red, F555W in green, and F336W in blue. Black-and- 
white images are the image cutouts in the F814W and the Ha 
filters of regions A-E, as outlined in green in the central panel. 
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F555W - F814W F555W - F814W 

Figure 14. CMD and color-color diagram of Wolf-Rayet sources in M83 marked in Figure 13. Three sources (78, 86, 105) circled in both diagrams are identified as 
clusters in our M83 WFC3 images. A twice solar metallicity BC09 cluster model (G. Bruzual & S. Chariot 2009, private communication; see also Bruzual & Chariot 
2003) is shown in the right panel (dots), and Padova stellar models (Marigo et al. 2008) are shown in both panels. No correction has been made for internal reddening, 
which appears to be minimal for all but one object. The arrow indicates the Galactic reddening vector. 



Figure 15. HST /ACS images (290" x 290", 68 pc x 68 pc) of the compact star cluster NGC 2108 in the LMC: a color composite of the F435W, F555W, and F814W 
ACS images (left), the F555W ACS image (middle), and the degraded F555W ACS image (right) by a factor of 100 (from a distance of ~50 kpc to a distance of 
~5 Mpc). See Appendix B for details. 


The original finding chart of region A taken from Hadfield et al. 
(2005) is shown in the top right panel of Figure 13. In the Ha 
image, the stellar continuum has been subtracted from the Ha 
image, leaving just the ionized gas. We matched the coordinates 
of the Wolf-Rayet sources and the HST images by assuming that 
the Wolf-Rayet source 99 is perfectly centered on a relatively 
isolated, compact Ha knot. We also compared the locations of 
the 17 matched Wolf-Rayet sources in the M83 HST images 
with the finding charts from Hadfield et al. (2005) to verify our 
source matching. 

The excellent spatial resolution of the HST image also makes 
it possible to check whether any of the Wolf-Rayet sources are 
actually compact star clusters. While the majority do appear to 
be individual stars, or a dominant star with a little “fuzz,” which 
is likely a faint companion star, there are at least three cases 
where the counterpart appears to be a bright compact cluster 
(i.e., 74, 86, 105; circled objects in Figure 14). These three 


objects are brighter than My — —9 mag with the concentration 
index (Cl) larger than 2.3, consistent with the cluster Cl used 
in Chandar et al. (2010). There are several more cases where 
the counterpart is in a looser association of stars (e.g., 78, 102, 
103). 

All the Wolf-Rayet sources, except the sources in cutout E 
and 108 in cutout D, are located in the spiral arms. This is 
expected, since these are the sites of most recent star formation. 
Nine of the 17 are in regions of strong Ha emission, seven more 
are in regions of faint Ha emission, and only one source (107 
in cutout D) appears to have no Ha associated with it. 

Our primary question is whether the Wolf-Rayet sources tend 
to be found in regions that we estimated to have young ages. In 
Table 2 we match the Wolf-Rayet sources with the region where 
they are found in and find that their ages range from 1 . 1 Myr 
(Region 35) to 50.1 Myr (Region 41). This shows that such a 
correlation does exist, with 6/8 (75%) of Wolf-Rayet sources 
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Figure 16. Image cutouts of HST /ACS F555W image of NGC 2108. The large circles in yellow show the 300 pixel aperture (corresponding to a ~3 pixel aperture 
at the distance of M83) centered at the location determined from the degraded image in Figure 15. The location of the dominant star (used to determine the “truth” 
value) is shown in red. See Appendix B for details. 


in regions having Age Reg <4.6 Myr. Similarly, Hadfield et al. 
(2005) found that five of their Wolf-Rayet sources are included 
in a compact star cluster catalog developed by Larsen (2004), 
and three of these five clusters (193, 179, 179) have estimated 
ages in the range 1.5-6 Myr, consistent with ages expected for 
Wolf-Rayet stars. 


Figure 14 shows the CMD and color-color diagram of the 
measured colors from our HST observations of the object that 
best matches the Wolf-Rayet sources from Hadfield et al. 
(2005). We find that most of these sources are in the top left of the 
two-color diagram. We conclude that there is a strong tendency 
for Wolf-Rayet sources to be associated with the youngest 
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regions of star formation. The good correlation provides general 
support for the accuracy of our CMD age estimates. 

APPENDIX B 

EFFECT OF SPATIAL RESOLUTION 

Spatial resolution affects virtually all studies of individual 
stars in nearby galaxies. As is typical, in this paper we have 
assumed that the point sources identified in Section 2 are 
individual stars, although we also discovered a small fraction 
of pairs of nearly aligned stars, one red and one blue, that 
can be identified from their anomalous colors (e.g., “pseudo- 
match” Figures 2 and 4). Here, we address more generally how 
spatial resolution affects the measured colors and magnitudes, 
and hence the derived ages, of stars in our M83 regions. 

We begin by answering a simple question: “What would the 
Orion Nebula look like at the distance of M83 ? Would we be able 
to distinguish the four central stars that make up the Trapezium, 
or would they appear as a single star?” The separation between 
the four stars in the Trapezium (10" at D 0.5 kpc) would be 
roughly O'.' 001 at the distance of M83, approximately a factor 
of 40 smaller than a single pixel in our image, and hence these 
would appear as one point-like source. However, we note that the 
vast majority of the remaining bright stars in the Orion Nebula 
are much more widely separated and would not suffer from this 
problem. We also note that our selection criteria are designed 
to avoid including stars in the crowded central regions of star 
clusters and associations, which also minimizes the problem. 

We perform a more quantitative analysis by degrading HST / 
ACS images of the compact star cluster NGC 2108 in the LMC 
(ID: 10595; PI: Goudfrooij; see Goudfrooij et al. 2011), which 
is at a distance of ~50 kpc, by a factor of 100 (i.e., equivalent 
with the distance of »5 Mpc), as roughly appropriate for M83. 
The degraded images, created from a combination of rebinning 
and Gaussian smoothing to approximately mimic the PSF of 
Advanced Camera for Surveys (ACS), are then run through 
our normal object-finding software described in Section 2. The 
numbers of objects found in the degraded NGC 2108 images are 
~20 in the F555W and F814W ACS images. Figure 15 shows 
the F555W ACS image of NGC 2108, along with the degraded 
images and a color image (i.e., a color composite of F435W, 
F555W, and F814W ACS images). 

Additionally, we then perform aperture photometry on the 
detected source positions in the non-degraded ACS images 
using a radius of 300 pixels to collect light from all stars that 
would fall within our aperture if NGC 2108 was located at 
a distance of 5 Mpc (i.e., M83). Objects along the edge and 
within 1000 pixels of the center of the cluster in NGC 2108 
are discarded, as would be the case for the corresponding 
photometry in M83. Figure 16 shows the extent of the 300 pixel 
aperture for several sources in yellow, as well as the nearest 
dominant star in red. The magnitude of the dominant stars (using 
a 3 pixel aperture and the appropriate aperture correction) is 
considered the “truth” measurement, while the magnitude in 
the corresponding 300 pixel aperture is the value that would be 
determined at the distance of M83. We calculate the magnitude 
and color offsets, i.e., A I and A( V — I) mag of the dominant 
stars, by taking the difference between the 3 pixel and 300 pixel 
aperture magnitudes. 

One complication is that our default sky subtraction method 
would remove the vast majority of faint and moderate brightness 
stars in the background annulus in the NGC 2108 image, 
but these stars would not be detected individually at the 



F555W - F814W 

Figure 17. Original CMD of Region 3 in M83 (solid points) with the effect of 
simulation representing a factor of 100 degradation in spatial resolution (open 
points). The Age CMDa uto of this region changes from 21.9 Myr (determined 
from the solid points) to 23.8 Myr (from the open points). The difference in the 
age estimates is within ~la error listed in Table 1. See Appendix B for details. 


distance of M83 and hence would be included in the sky 
measurement. We therefore adjusted the “clipping” parameters 
in the sky subtraction algorithm to remove only the bright stars 
in NGC 2108, thereby mimicking the case for M83 to the degree 
possible. 

In Figure 17 we show the original CMD for Region 3 in M83 
(solid points), compared with a version where the magnitudes 
and colors have been “perturbed” by an amount A I and A (V — I) 
determined by matching to the closest values of (V — /) from 
the NGC 2108 data experiment (open points). As expected, 
corrections for the degraded spatial resolution in M83 tend 
to make the corrected magnitudes slightly fainter and redder, 
as can be seen in Figure 17. The mean difference in V is 
—0.24 ± 0.30 mag with a range from +0.20 to —0.65 mag. 
In /, the mean difference is —0.40 ± 0.39 mag with a range 
from +0.07 to —1.18 mag. We note that the largest difference 
is seen for object 10, as shown in Figure 16, where several 
relatively bright red stars in the 300 pixel aperture cause the 
largest correction (i.e., —1.18 mag in I and redward by 0.63 mag 
in (V — /)). Object 40, dominated by a single very bright star, 
shows the opposite extreme, with a change of —0.02 mag in 
V and —0.01 mag in I. We also note that object 31 does not 
make it into our sample since the two bright stars are far enough 
apart that it is clear from the degraded image in Figure 15 that 
this is not a single star. This object gets removed in the DoPHOT 
photometry, as described in Section 2, and hence is not used in 
the experiment described above. 

We then rerun our age-dating software on the corrected CMD 
in Figure 17. We find a mean age of 23.8 Myr compared to the 
original value of 21.9 Myr. Hence, while a specific “star” may 
be affected by a sizable amount, the overall effects of degraded 
spatial resolution are relatively minor for our study. 
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To summarize, we have performed a numerical experiment 
using observations of an intermediate-age star cluster in the 
LMC (NGC 2108) for our “truth” image. Using a factor of 100 in 
spatial degradation, roughly appropriate for our measurements 
in M83, we find that the typical corrections to our photometry 
are on the order of a few tenths of a magnitude, although much 
larger values are possible in specific cases. We conclude that 
while spatial resolution can result in measurable differences in 
the luminosities and colors of bright stars at the distance of M83, 
it does not strongly affect the age estimates for the 50 regions 
studied in this paper. 
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